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Abstract 

' We present a lattice calculation of the interquark potential between static quarks in a "full" QCD 

_ , simulation with 2 flavours of dynamical Wilson-quarks at three intermediate sea-quark masses. We 

work at f3 — 5.6 on lattice size of 16 3 x 32 with 100 configurations per sea-quark mass. We compare 
the full QCD potential with its quenched counterpart at equal lattice spacing, a -1 ~ 2.0 GeV, which 
is at the onset of the quenched scaling regime. We find that the full QCD potential lies consistently 
below that of quenched QCD. We see no evidence for string-breaking effects on these lattice volumes, 
V ~ (1.5 fm) 3 . 

H— > . 

i ! 1 Introduction 

Ph. 

Ever since the seminal paper of Creutz|IJ on the confining character of the static quark-antiquark in- 
teraction in the weak coupling regime of quenched QCD, refined lattice methods have been devised and 
, ' targeted towards an improved determination of the interquark potential. By now, quenched studies on 

large lattices (> 32 4 ) have reached a statistical accuracy in the percent region at small and medium 
interquark separations, i?, for bare lattice couplings in the scaling regime 6.0 < (3 < 6.8@, §, f§. This 
corresponds to lattice resolutions in the range 1.94 < a -1 < 6 GeV, which is fine enough to observe 
running coupling effects (jjJQ. 

As full QCD lattice simulations with Wilson fermions are presently still in their infancy it is rather 
obvious to ask the question whether the form of the static potential can provide a decent signal for 
fermion loop effects in the real QCD vacuum state. The shape of the static potential is expected to 
grow more convex as effects from dynamic fermions turn gradually stronger with decreasing quark mass; 
however, this can only be observed provided we can access a large enough window of observation in R at 
sufficiently small quark masses. We can estimate the effect of dynamical fermions on the Coulomb-term 
of the potential from the lowest-order coupling constant : 



aiR) ^-tb \ogRA v ' (1) 
where the iVy-dependence is in bo (Nf denotes the number of flavours): 

b °- 48^ ■ ( 2 ) 

1 At present, such level of precision is extremely hard torcalizedn lattice evaluations of hadronic attributes proper, even 
in quenched QCD, since this requires very high statisticsls, H, 
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Naively, from bo only, we would expect to observe an effect on the order of +14 % when switching from 
Nf — to Nf = 2. Obviously, to be sensitive to unquenching, one requires data with an overall error 
(statistical and systematic) below the 5 % level. Such a precision was clearly not reached by previous 
lattice studies with dynamical Wilson fermions [ll| in the scaling regime, which we expect to start 

at (3 ~ 5.6.0 

In this paper, as part of our ongoing 100 Teraflopshrs simulation of 2-flavour Wilson-fermion QCD, 
which was described in Ref . []l3f , we will present a first high statistics study of the interquark potential. 
We work at a bare lattice coupling of (3 — 5.6 on a lattice volume of 16 3 x 32 and simulate at three 
sea-quark masses which correspond to the hopping parameter values K sca — 0.156,0.157 and 0.1575. 
With these parameters the lattice resolution is expected to be around a -1 ~ 2 GeV jl(J |l3|, 0. 

Field configurations are produced on Quadrics QH2 machines using the Hybrid Monte Carlo Algo- 
rithm; for details we refer the reader to a forthcoming publication(ll|. In this analysis, we exploit 2500 
trajectories per sea-quark massQ, which corresponds to three times 100 independent configurations (see 
section ||) . 

2 The Static Potential 

Autocorrelation The statistical quality of a sample of configurations generated as an HMC timeseries 
is largely determined by the integrated autocorrelation time T in t. We have obtained clean signals both for 
exponential and integrated autocorrelation times of the plaquette using the complete sets of trajectories. 
In table | we list the values for the integrated autocorrelation time Tint- Since the plaquettes are con- 



n sca 0.156 0.1570 0.1575 
n n t 4.8(7) 5.1(8) 14.8(8) 

Table 1: Exponential and integrated autocorrelation times of the plaquette. 



sidered to represent the worst case with respect to autocorrelations we choose to perform measurements 
on configurations separated by 25 trajectories (a detailed discussion of autocorrelation issues follows in 
a forthcoming publication)^]. 

Analysis Method The static (spin independent) potential V(R) is computed in a standard fashion 
from the path-ordered products of link variables around space-time rectangles. In order to enhance 
the ground-state signal, we use a local gauge-invariant APE-type fl(| smearing procedure on the spatial 
links. The smearing parameter a is set to a = 2. We have optimised the number N of smearing iterations 
and find N = 20 or N = 25 to yield overlaps > 80 % for all R with very little dependence on R (lattice 
quantities are hatted to avoid confusion). 

Similarly to the quenched case 0] , off-axis measurements are included for improved spatial resolution; 
this also allows an investigation of rotational invariance restoration. We perform measurements at 36 
different values of R comprising 8 on-axis and 24 off-axis constellations. 

The potential is defined through local masses at time T : V T {R) = log y^J^f ^ ■ We search for 

2 In this situation, one is tempted to assess the impact of dynamical fermions by matching, at given values of /3 and 
K sea , the Wilson lo ops W(R X T) onto their respective quenched analogues, by suitable /3-shifts, Af3(R X T, /3, K se a)- Using 
this technique ref. [hcj indeed found an increasing A/3 with R. However, as this procedure can only be related to physical 
quantities in the large T limit one would prefer to investigate the QCD potential directly Jl^ ] . 

3 This is approximately half of our scheduled sample. 

4 Note that due to a variable step size (A^J 111 = 100 and Af™J = 20) these configurations are not equidistant in 
molecular dynamics time. 
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plateaus defined by 



(3) 



V T {R) 

Good plateaus are found for T > 3. We have checked that local masses, V T (R), and semi-local masses, 
from constant fits of V T on three time-slices T = 3, 4, 5, yield the same values for the potential. Fur- 
thermore, the central values of the fit-parameters (as given below in tables |2| and ||) remain virtually 
unaffected: however, using semi-local masses one may reduce the errors by approximately 15 % (albeit 
at the cost of a slightly increased x 2 per degree of freedom). 

Lattice artefacts in the potential due to the difference of the one-gluon exchange on the lattice and 
in the continuum are corrected for by subtracting a term proportional to 

SV(R) = ([h I) , (4) 



R R 

where [i] denotes the lattice one gluon exchange extrapolated to infinite volume [Q, [l?], |l8| ]: 

[±] =A-KYim v ^G{R;V) with G(R;V) = Jk^fL > ( 5 ) 

where h = ^-rrii, m, = -L/2 + 1, . . . , L/2, V = L 3 . 

In order to fix the size of the correction term g SV we interpolate our data using the following four- 
parameter ansatz (the detailed form of the fit is of no interest here, we only wish to show that a smooth 
interpolation can be obtained): 

V(R) = V + kR- 4 +gSV(R) . (6) 
R 




Figure 1: (V(R) — gSV(R)j as a function of the lattice separation R. 

Results The data and fits are shown in figure [l] where we plot (v(R) — gSV(R) \ as a function of the 
lattice separation R. Throughout the analysis we exclude the two smallest R values from the fit^|. The 

5 These points are, however, included in figure EI 
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Figure 2: Relative deviations (V(R) — Vf(R))/Vf(R) of the potential values from the fit curves Vf(R), 



eq. 



resulting parameters are given in table ^|; errors are obtained from a bootstrap sample of size 250 and 
correspond to 68 % confidence level. We find all \ 2 P er degree of freedom to be around 1. 



^sca 


^0 




k 


e 


9 


X 2 / d.o.f 


0.156 


0.724; 


f r 

- 6 


0.0514±^ 


0.324 j; I 


0.361 2 


18/30 


0.157 


0.721; 


f 8 
- 6 


0.04451}? 


0.3161^ 


0.351 2 


31/30 


0.1575 


0.728; 


f 7 
- 7 


0.03961}? 


0.323±ig 


0.371 2 


25/30 


quenched 




6.0 


0.6591 


+ 43 
-43 


0.04771 


0.29171^ 


0.29161^4 


10/29 


6.2 


0.6336 


+ 24 
-24 


0.02574±^ 


0.2865 tit 


0.2651}} 


43/63 



Table 2: Fit results from the ansatz eq. |6[ For a discussion of the quenched data see section | 



A more detailed view of the quality of the data and the parametrisation can be obtained by looking 
at the relative deviations of the data points from the interpolating curves as shown in figure We find 
the deviations to be smaller than ±2% with minor exceptions at very large R. We conclude that adding 
a term proportional to 5V(R) is sufficient to remove lattice artefacts at small R and that our interquark 
potentials show no sign of non-linear behaviour towards large R, as predicted from string breaking. 



3 Scaling the data 

In the case of full QCD it is appropriate to find the lattice scale from a quantity defined at an intermediate 
distance rather than using the asymptotic quantity "string-tension". Therefore, following ref. p|, we 
extract the lattice spacing ap by matching the force from our lattice potential to the force obtained 
from potential models as applied to quarkonium spectroscopy. This amounts to finding the numerical 
solution, in terms of Ro, of the equation 

R 2 o% = 1-65 , (7) 

ClR \R„ 
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which, in phcnomcnological models, is given by the reference length 



Ro = a F R Q = 0.5 fin . (8) 

Since this procedure involves the use of a numerical derivative it is clearly more delicate than setting 
the scale by the string-tension and hence requires precise data. 
The force values are calculated from the finite differences 

f(r) = , (9) 

' ' (Ri - R 2 ) 

where V C (R) are the corrected data which we view as our lattice estimates of the continuum QCD 
potential : 

V C (R) = V(R) - gSV{R) . (10) 

All possible pairs of potential data with distances 1.4 > (i?i — R2) > 1 are used to calculate an improved 
R in eq. |9| by demanding that 

-^^=F(R), (11) 

where 

V theo {R) = V + kR-^ . (12) 

Eq. [ll] is then solved for R using the fit parameters of table 0. 

The resulting values Rq and lattice spacings a? are given in table ||. We find very little difference be- 
tween the values of Rq obtained from this numerical difference procedure and those obtained analytically, 
with the fit parameters of eq. || : 

This is taken as further evidence of the smoothness of our data and as an indication for the stability of 
the method. Note that the statistical errors are as low as 3 %. 



K>sc& 


Ro 


a^[GeV] 




0.156 


5.131 I 


2.02 1 I 


1.941 I 


0.157 


5.55 -yj 


2.191 \ 


2.091 % 


0.1575 


5.83^2 


2.30 ± 5 


2.2ll I 


quenched 




6.0 


5.35 1 3 


2.111 \ 


2.02 1 I 


6.2 


7.27± I 


2.861 i 


2.741 I 



Table 3: Lattice spacings obtained from the force at Rq = 0.5 fm and from the string tension k = 
a 2 x 0.1936 GeV 2 (see table | for k). Err ors are statistical only. The quenched data are discussed in 
section 0. 



As we see no deviation of the potential from a linear behaviour for large R we can also attempt to 
extract the lattice spacing from the force at larger distances, i.e. at around 1 fm, where our data are 
dominated by the linear term in eq. [(]. In practical terms, we then interpret the values of k from table || 

as string-tension values from which we derive the lattice spacings ay = \J k/a, where we take, somewhat 
arbitrarily, ^fa = 0.44 GeV (see table p|). Note that the sizes of the statistical errors are beginning to 
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be competitive to those of previous quenched calculations (without link- integration) j|, |2Cfl. Moreover, 
we observe that the ratio of/ov remains unchanged for all three sea-quark masses indicating that a 
slightly different choice of y/a, which we find to be 458 ±12 MeV for n = 0.1575, renders consistent 
lattice spacings from both schemes. 

We comment that our hadron spectrum analysis yields lattice spacings consistent with those of table 
||, when using the mass of the p to set the scale Q . 

In figure || we show our results in scaled form. We find that the data collapse to a universal potential 
over the complete range of our measurements, R € [0.1, 1.5] fm. An analysis similar to that of figure ^ 




0.5 1.0 1.5 

R[fm] 



Figure 3: Scaled potential (using ay) in physical units, 
shows that nearly all points lie within ±2 % of the universal curve. 

4 Discussion : effects of dynamical fermions ? 

At first sight, it might seem disappointing to find that our scaled data collapse to a universal function; 
clearly, this shows that the variation of K sea from 0.156 to 0.1575 is too small to resolve mass-dependence 
in the static potential. However, we can do more/better by comparing our data to the extreme case, the 
case of infinitely heavy sea quarks, that is, the quenched approximation. 

Table [3] shows that the lattice spacings in the present calculation correspond to quenched couplings 
slightly above (3 = 6.0. We therefore performed two high statistics runs in pure gauge theory including 
link-integration at j3 = 6.0 (lattice size 16 4 , 570 configurations) and (3 = 6.2 (lattice size 32 4 , 100 
configurations) ||T), the j3 — 6.2 simulation being done to investigate scaling. In order to achieve 
comparable systematic errors the evaluation of the quenched potential data and their analysis follow the 
same procedure as described in section 2 and we perform four-parameter fits at both /3-values using a 
common lower cut R > 0.196 fm, designed to give good x 2 /d.o./. for all fitsQ. 

6 We have checked however, that varying the fit range to include more values of V(R) at smaller R leaves our conclusions 
largely unaffected. The fit parameters for the full QCD case are, in fact, completely stable. 
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Results are shown in tables |2j and |3|. We find the quenched data, whose lattice spacings vary by 
about 40 %, to fall upon a universal curve with deviations smaller than 1 %, corroborating the earlier 
results of ref j|. Furthermore, since the lattice spacings at /3jv,=o = 6.0 roughly correspond to those of 
PN f =2 — 5.6 (see table ||) we ar e indeed at the onset of the scaling regime for our full QCD simulation. 

The comparison of quenched and unquenched QCD is shown in figure ^ where we plot the data in 
the small i?-region, R < 0.4 fm. We find the full QCD data to lie significantly below that of quenched 
QCD ! In this scaling representation the data are solely sensitive to the strength of the Coulombic term 
e. In order to quantify the effect of dynamical fcrmions we can therefore directly extract the information 
from the e-values as quoted in table |[ As for both cases, full and quenched QCD, we observe good 
scaling, we can represent each data-set by its appropriate average value, ejv /= 2 = 0.3210(100) and 
ejv.=o = 0.2890(55). From these numbers, we find unquenching to induce a decrease of 11 % in the 
coulombic part of the static potential. 




Figure 4: Comparison of scaled potentials from quenched and full QCD (combined data for quenched 
((3 — 6.2 and (3 = 6.0) and full (three sea-quark values) QCD). The plot shows both fits with their 
one-sigma errorbands where only points right of the vertical line were used for fits, R £ [0.196, oo] fm. 
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5 Summary and conclusions 



We have presented a calculation of the static potential which is precise enough to unravel, for the first 
time, the effects of unquenching at intermediate distances. We find the full QCD coupling constant to lie 
consistenly above that of quenched QCD. However, the data are not sensitive enough to see variation of 
the potential with the sea-quark mass. Our method is based on a consistent modelling of the Nf — 2 and 
Nf = potentials. The full QCD lattice spacings are in the range a -1 ~ 2 GeV which, in the quenched 
case, is at the onset of the scaling regime, /3jv /= o — 6.0. Our result is in accord with the JV/-dependence 
of a as extracted from the perturbative expansion of the plaquette (see for example |23] ) . 

The lattice volume of (1.5 fm) 3 was not large enough to observe string breaking. Note that using 
quenched data, one can obtain an upper bound for the distance at which string breaking should be 
observable to be at 1.7 fm [M. 

We are currently performing a simulation with the same lattice-resolution but 50 % extended lattice 
size, in order to assess finite-size effects and to work at lighter quark masses [^2| . 

We are encouraged by these results in our systematic search for Sea quark Effects on Spectrum And 
Matrix elements. 
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